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ABSTRACT 

Extremely high deuteration of several molecules have been observed around low mass protostars 
since a decade. Among them, formaldehyde and methanol present particularly high deuteration, 
with observations of abundant doubly and triply deuterated forms. Both species are thought to be 
mainly formed on interstellar grains during the low temperature and dense pre-collapse phase by 
H and D atom additions on the iced CO. We present here a theoretical study of the formaldehyde 
and methanol deuteration obtained with our gas- grain model, GRAINOBLE. This model takes into 
account the multilayer nature of the mantle and explores the robustness of the results against the 
uncertainties of poorly constrained chemical and surface model parameters. The comparison of the 
model predictions with the observations leads to two major results: i) the observed high deuteration 
is obtained during the last phase of the pre-collapse stage, when the density reaches '--^ 5 x 10^ cm~^, 
and this phase is fast, lasting only several thousands years, ii) D and H abstraction and substitution 
reactions are crucial in making up the observed deuteration ratios; This work shows the power of 
chemical composition as a tool to reconstruct the past history of protostars. 



1. INTRODUCTION 

Although the deuterium e lemental abun dance is 1.5 x 
10~^ relative to hydrogen (jLinskvl l2003f ). observations 
carried out during the last decade have revealed high 
abundances of singly, doubly, and even triply deuter- 
ated molecules in low mass pre-stell ar cores and Class 
protostars (jCeccarelli et al.l 120071 ). In this work, 
we focus on formaldehyde and metha nol, both species 
likely synthesized on the grain su rfaces (jTielens fc Hagenl 
fl98a IWatenabe fc Kouchil[200l as gas phase reactions 
are unable to reproduce the large obse rved abundances 
(|Roberts et al.|[20p : lMaret et alj[2005l ). 

In pre-stellar cores, HDCO/H2CO and D2CO/H2CO 
abundance rati os up to 10% have been observed 
(jBacmann et al.l 120031 : iBergman et al.l 120111 ). The ob- 
served deuterium fractionations increase with the in- 
creasing CO depletion, suggesting that the latter is 
a key parameter. Analogously, several formaldehyde 
and methanol isotopologu es have been detected i n Class 
low mass p rotostars (iCeccarelli et al.l 119981 120011 : 
iParise et al.ll2002 . _200i, |2006D. However, unlike in pre- 
stellar cores, no correlation between the CO depletion 
and deuterium fractionation is observed, so that the 
deuteration process is thought to occur in the previ- 
ous pre-collapse phase. In general, methanol is more 
enriched in deuterium than formaldehyde. In most of 
the observed sources, the simply deuterated molecule 
HDCO shows a deuterium fractionation between 13 and 
20% whereas CH2DOH shows a fractionation between 
37 and 65% (IParise et al.ll2006l) . CHD2OH also shows 
a higher fractionation than D2CO, by a factor of ^ 2. 
Finally, the deuteration process occurs more efficiently 
on the methyl group than on the hydroxyl group of 
methanol, as the [CH2DOH]/[CH30D] ratio has been ob- 
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served to be between 10 and 20 in a sample of low- and 
intermediate- mass pro tostars bv lParise et al.l ()2006[ ) and 
iRataiczak et al.l ()201lD . 

Recent experimental works have confirmed the syn- 
thesis of formaldehyde and methanol by hydrogena- 
tion of iced CO (Watana be & Kouchi2002,: Hidaka et all 
120071: iFuchs et al. 2009) and also highlighted the com- 
plex chemical evolution of H2O-CO ices when th ey are 
expos ed to D and H atoms. In deed, [N agaok a et al.l 
(|2005f ). iNagaoka etall (|2007| ) and iHidaka et al.i (|2009l ) 
have shown that formaldehyde and methanol can be 
efficiently deuterated into HDCO, D2CO, CH2DOH, 
CHD2OH and CD3OH, when they are irradiated by 
D atoms. Conversely, abstractions of D can only oc- 
cur on formaldehyde, and not on methanol. Therefore, 
formaldehyde isotopologues only can be hydrogenated to 
form back HDCO and H2CO if they are exposed to H 
atoms. The relative reaction rates deduced from these 
works show that H and D abstractions and substitutions 
on formaldehyde and methanol are as efficient as addi- 
tion reactions ( see also the theoretic al calculations by 
IGoumansll20TTI: IGoumans et al.ll20TTI ). These processes 
could, therefore, largely increase the deuterium fraction- 
ations of formaldehyde and methanol after the complete 
depletion of CO and their formation on the surfaces. 

In the past, astrochemical models have struggled to re- 
produce the observed deuteration ratios of formaldehyde 
and methanol. While it is now clear that CO depletion 
plays a major role in increasing the ato mic D/H ratio 
of th e gas landing on the grain surfaces (|Roberts et al.l 
I2004D , a full model coupling the gas and grain chemistry 
that simultaneously reproduces the observed formalde- 
hyde and methanol deuteration is still missing. A p revi- 
ous attempt was carried out bv lCaselli et al.l ()2002t ) and 
Stantcheva & Hcrbst (2003), who studied the formalde- 
hyde and methanol deuteration as function of the atomic 
D/H ratio, ta ken as a free param eter. However, as also 
emphasised bv lParise et al.l ()2006[) . the use of a constant 
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[D]/[H] ratio is i) unable to predict all the deuterium 
fractionations at the same time and ii) is not necessarily 
correct, as formaldehyde and metha nol may b e formed on 
the grains at different times. Caza ux et al.l ([2011,) have 
studied the deuteration of water and formaldehyde (and 
not methanol) by coupling gas phase and grain mantle 
chemical networks. 

In this article, we re-consider the problem of the 
formaldehyde and methanol deuteration, using our gas- 
grain coupled model, GRAINOBLE, that takes into ac- 
count D and H atoms addition, but also abstraction and 
substitution reactions. Our goal is to simultaneously re- 
produce the formaldehyde and methanol deuteration ob- 
served towards Class protostars by considering their 
formation on interstellar grains only. Note that, once 
sublimated from ices, the abundance of these species is 
only slightly affected by gas phase reactions, as the typ- 
ical chemical timescale (~ lO'^ yr, ICharnlev et al.|[l997[ ) 
is l a,rger than the typ ical age of these objects (~ 10'' 
yr, I Andre et al.ll2000l ). We will show that understand- 
ing how and when the observed deuteration occurs will 
also provide us with hints on the process itself and on 
the past history of the protostars. 

2. MULTILAYER MODELING OF DEUTERATED ICES 

2.1. Description of the GRAINOBLE model 

For this study we use the GRAINOBLE model 
(|Taauet et all 1201 IL hereafter TCKll). Briefly, it is a 
gas-grain chemical mo del, based on the rate equations 
approach introduced bv lHasegawa et afl ()1992f ) for time- 
dependent grain surface chemistry modeling. It considers 
the following four processes: 

i) Accretion of gas phase species onto the grain surfaces 
as function of time. 

ii) Thermal diffusion of adsorbed species. The hopping 
rate follows a Boltzmann function which depends on the 
diffusion energy Ed- 

iii) Surface reactions via the Langmuir-Hinshelwood 
mechanism. The probability of the reaction is given by 
the exponential portion of the quantum mechanical prob- 
ability for tunneling through a square barrier and there- 
fore depends on the activation energy Ea of the reaction. 

iv) Thermal evaporation caused by the thermal balance 
and by the cosmic-ray ind uced heating of the grains 
(jHasegawa fc Herbs^ll993al) . We ignored the photolytic 
process due to the cosmic-ray induced UV photons be- 
cause they have a negligible influence on the abundance 
of formalde hyde and methanol in dark clouds (5% at 
maximum, iCuppen et al|[200l . 

We follow the mantle formation on grains with a mul- 
tilayer approach in which the outermost layer only is 
reactive, while the mantle bulk remains inert. Unlike 
IHaseeawa fc Herbstl (I1993bl) and lGarrod fc Paulvl (l20TTh 
models, the trapping of particles into the bulk is per- 
formed one layer at a time, once the considered layer is 
filled. 

2.2. Chemical network 

We consider the accretion of gaseous H, D, H2, HD, 
D2, O, and CO onto grains. 

We assume that formaldehyde, methanol and their 
deuterated isotopologues are formed via hydrogenation 
and deuteration addition reactions on molecules initiated 



by the accretio n of CO. Based on experimental and the- 
oretical works (iHidaka et al.ll2007t iWatanabe fc Kouchil 
[2008t lA^rsson et al.ll2'oTT[ ). the energy barriers of the 
D and H addition reactions involving CO and formalde- 
hyde are very similar, so we assume them to be identical 
(= Ea). Note that, however, given the high uncertainty 
in the value of Ea (see TCKll), this is taken as a free 
parameter. In contrast, reactions involving a radical are 
barrier less. 

Following the experimental works of IHidaka et al.l 
(2009) and Nagaoka et al. (2007), we also include the 
abstraction and substitution reactions of H and D on 
formalde hyde and methanol accordin g the scheme pro- 
posed by IWatanabe fc Kouc"hil ()2008f) and IHidaka et al.l 
(|2009f) . and shown in Fig. [TJ We adopt the probability 
of each reaction following the relative rates deduced ex- 
perimentally when they are available, as marked in the 
Figure. For the reactions that have not been derived by 
the experiments, we adopt the probability measured for 
the same isotopologue by analogy. For example, the re- 
action D2CO -I- D — CD3O has a reaction rate equal to 
the H + CO rate muhiplied by 0.1 (D + CO) and 0.66 
(D2CO + H). 

Even if most of water is likely formed during the 
translucent phase, and therefore before H2CO and 
CH3OH fsee lOberg et al.l 1201 ID . our model also takes 
into account its formation, as it is in competition 
with formaldehyde and methanol formation according 
to TCKll. With respect to TCKll, we added the 
deuteration reactions and the path 

H2 + OH ^ H2O + H 

which seems to be the most efficient reaction in m olecular 
clouds, as suggested bv lCuppen fc HerbstI (|2007D . 

Since the deuteration of formaldehyde and methanol 
on the grains depends on the atomic gas phase D/H ra- 
tio, we are here particularly interested in the chemical 
network leading to the formation of atomic deuterium in 
the gas p hase. We use the f ractio n ation reactions intro - 
duced bv lRoberts fc Millaii (f2000l) : iRoberts et al.l pOOl 
to study the formation of Hj^ isotopologues and atomic 
deuterium. The gas phas e chemical ne t work is shown 
in Figure [T] As shown by iFlower et al.l ()2006[ ) , the or- 
tho/para ratio (opr) of H2 can influence the deuteration 
of H^ and consequently the atomic D/H ratio, but only 
when it is > 10"'^. The available measures in cold gas in- 
dicat e low H2 opr values (< 10""^ — 10~^, IDislaire et al.l 
120121 and references therein). We, therefore, did not 
consider this effect in this work. 

2.3. The physical model 

Rather than simulating the evolution from the dif- 
fuse cloud state to the pre-stellar core, we focus on the 
last stage of the evolution, when the material is already 
molecular and the density reaches the value nu ■ The ini- 
tial abundance of gas phase species (H, D, HD, D2, H^^ 
isotopologues, ...) are given by the steady-state abun- 
dance, obtained solving the gas phase chemical network 
presented in §2.2 and considering the recombination of H 
and D onto grains. The abundance relative to H nuclei 
of the deuterium and CO reservoirs are 1.5 x 10~^ and 
4.75 X 10~^ 



3 




53 



{ HCO* ^^^^ f^ vH^^ H+H, ) 



CD,0,CD,OD CH,OH,CH,0 
' ' ^ CHD0D,CD20H CH^OD, CHDOH ^ ' ' 




DCO* 



+H2 
+C0 



H,D< 



+HD 

n +e, gr 



+ H;,H 

+C0 r 



+HD,D 

~1 gr 



+Hj,H 
+C0 r 



+HD,D 
n +e, gr 



CD3OD ICD3OHI CHD2OD [cHD;OH| CH;DOD rai;DOH I icHjODl ICHjOhj 



Fig. 1. — left) chemical network on grain mantles focusing on the formation of formaldehyde and methanol. Solid, dashed, and open 
arrows refer to addition, abstraction, and substitution reactions. Black arrows refer to reactions involving H while grey ar rows represent 
reactions involving D. Finally, the bold values refer to the rates relative to the H + CO reaction experimentally deduced bv lNagaoka et all 
((2007) and Hidaka et al. (2009) while non-bold values show the deduced relative rates by analogy with similar reactions, right) Chemical 
network in the gas phase leading to atomic deuterium formation. 

We follow, then, the formation and evolution of the 
grain mantles during this phase, keeping the gas density 
hh constant. Briefly, as the time passes, CO freezes 
out onto the grain mantles, where it forms H2CO and 
CH3OH and their deuterated isotopologues. Therefore, 
the deuteration of the formaldehyde and methanol on the 
mantle primarily depends on the gaseous atomic D/H 
ratio. 



2.4. Multi-parameter approach 

We use a multi-parameter approach as described in 
TCKll, to take into account the variation or the uncer- 
tainty of physical conditions and surface parameters. 

Low mass pre-stellar cores are the place where most of 
formaldehyde and methanol are assumed to be formed. 
These objects show spatial distributions of densities of 
H nuclei between 10'' and 5 x 10^ cm~^ an d tempera- 
tures between 8 and 12 K (see lCrapsi et al.|[2007 1. Ac- 
cordingly, we consider here four density values (nn = 
10'', 10^, 10^, 5 X 10^ cm~^) and three temperatures 
(Tg — Td — 8,10,12 K). Besides, prestellar cores show 
an increase of the grain sizes caused by the coagulation 
process whose the efficiency increases with the density 
(|Flower et al.ll2005t fVastel et al.ll2006( ). We adopt three 
grain sizes aa = 0.1,0.2,0.3 /xm, which affect the deple- 
tion rates and therefore the density of H and D in the 
gas phase. 

As discussed in TCKll, the diffusion to desorption en- 
ergy Ed/Ef, ratio mainly depends on the ice properties. 
We, therefore, vary Ed/Et between 0.5 and 0.8, as sug- 
gested by experimental studies. The activation energy 
Ea of the reactions involving CO and formaldehyde is 
also a free parameter. By comparing the model predic- 
tions with the observations of solid CO and methanol, we 
have deduced in TCKll that Ea must be lower than 1500 
K, so that we restrain the value of Ea between 400 and 
1400 K. Light particles (H, D, H2, HD, D2) that accrete 
onto the ices have a binding energy distribution that de- 
pends on the ice properties, and the adsorption condi- 
tions of the adsorbate (see for example iHornekasr et al.l 



l2005f l. We therefore consider the binding energies rela- 
tive to amorphous water ice of these light species as a 
free parameter whose values are: 400, 500, 600 K. 

As the abundance of the gaseous atomic oxygen in the 
pre-stellar cores is uncertain, we consider three values of 
oxygen abundance relative to H nuclei X{0)mi- 10~^, 
10-^ 10-4. 

Finally, in TCKll, we have shown that the site size dg 
and the porosity factor Fpor have a little impact on the 
chemical composition of the grain mantles, so we assume 
ds = 3.1 A (iJenniskens et al.lll995[ ). and Fpor — 0. 

We run a grid of 2916 models in which we vary the 
parameters described above and listed in Table [TJ For 
each density nn, we compute the mean fractionation of 
each isotopolo gue with i t s 1 si gma standard deviation, 
following Wakelam et al.l (|2010l) . 



TABLE 1 

List of parameters and the values range explored in this 

WORK. 



Parameter 



Values 



Density nu 


10" - 10^ - 


lO'^ - 


5 X 10** cm-^ 


Temperature Tg = 


8 


10 - 


12 K 


Grain size 


0.1 - 


0.2 - 


0.3 /im 


Energy ratio E^/E), 


0.5 


- 0.65 - 0.8 


Binding energy Et^^atCH) 


400 - 


500 


- 600 K 


Activation energy Ea 


400 - 


900 


1400 K 


Initial oxygen abundance X{0)i„i 


10-8 


- 10- 


6 - 10-1 



3. RESULTS 

3.1. The D/H gas phase ratio 

At first approximation, the steady state densities of H 
and D, assumed as initial conditions, are roughly con- 
stant regardless the total density and their abundance 
relative to H nuclei, therefore, decrease with increasing 
riH- At low densities, a significant fraction (~30%) of 
the deuterium reservoir is already in atomic form before 
the depletion of CO: the increase of the atomic D/H ra- 
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tio with the CO depletion will therefore be low. On the 
contrary, at high densities, only a negligible fraction of 
deuterium is in the atomic form at the beginning, so that 
the atomic D/H ratio strongly increases with the CO de- 
pletion. The larger the density, the larger the gaseous 
atomic D/H ratio increase with the CO depletion, as 
highlighted in Figure [5] which shows the evolution of the 
atomic D/H as function of CO depletion, for tih equal 
to lO'^ and 5 x 10^ cm respectively. 

3.2. Model with addition reactions only 

In this section, only addition reactions are considered 
on grain mantles (see §2.2). In this case, the deuteration 
of formaldehyde and methanol primarily depends on two 
factors, as illustrated in Fig. [2| i) The increase of the gas 
phase atomic D/H ratio with the CO depletion: the ra- 
tio increases with increasing CO depletion, namely with 
time, and with increasing density (see above); ii) When 
formaldehyde and methanol are formed: as explained in 
detail in TCKll, the increase of nn delays the forma- 
tion of formaldehyde and methanol. In general, the two 
effects results in a larger deuteration of formaldehyde 
and methanol for larger densities and larger evolution- 
ary times. 




fp(CO) 



1000 



Fig. 2.— Normalized H2CO (solid lines) and CH3OH (dashed 
lines) formation rates (top panel) and gaseous atomic D/H ratio 
(bottom panel) as function of the CO depletion factor foiCO) = 
ng{CO)/ng^ini(CO), for njj = 10^ (thin red lines), and 5 X 10® 
(thick green lines) cm"''. The spikes are caused by the multilayer 
nature of the mantle. 



This is illustrated in Fig. [3J which shows the tem- 
poral evolution of the mean deuterium fractionations 
of iced formaldehyde and methanol along with their 
uncertainty, for different densities. The high density 
cases {uh > 10^ cm^'^) show large enhancements of the 
CH2DOH/CH3OH and HDCO/H2CO ratios with time, 
reaching the unity for the singly deuterated forms. No 
increase of the deuterium fractionation is, on the con- 
trary, observed at low densities {tih = 10^ — 10^ cm"'^), 
the deuteration ratios never exceeding 0.02. On the 
same figure, we re port the range of observed values (from 
IParise et al.l 120061 ). No model can reproduce the full set 
of observations. Indeed, the high density cases predict 
CH2DOH/CH3OH and HDCO/H2CO ratios in agree- 
ment with observations for times between 5 x 10^ and 
10^ yr. However, they cannot predict the observed high 
abundances of the doubly deuterated molecules. In ad- 
dition, observations show that methanol is three times 
more deuterated than formaldehyde whereas the models 
of Fig. [3]predict a factor 1.5 at most. We conclude that 
this class of models fails to reproduce the observations. 

3.3. Abstraction and substitution reactions 

Figure 0] shows the temporal evolutions of formalde- 
hyde and methanol deuterium fractionations when ab- 
straction and substitution reactions are included. From 
Fig. m it is clear that the inclusion of the addition and 
substitution reactions strongly increase the deuterium 
fractionations, especially at high densities. For densi- 
ties '--^ 5 X 10^ cm~^, the model predicts fractionations 
of doubly-deuterated molecules compatible with obser- 
vations, at time ~ 5 x lO'^ yr. Unlike CH2DOH, which 
shows an increase of deuteration compared to Fig. [3] 
and more particularly at longer timescales, the HDCO 
deuteration is not enhanced. This is due to the effi- 
ciency of deuterium abstraction on HDCO which allows 
the formation of H2CO whereas deuterium abstraction 
reactions on CH2DOH, leading to CH3OH formation, 
are negligible. The observed [CH2D0H]/[HDC0] ratio 
of 3 can now be predicted. Finally, because the abstrac- 
tion reactions can only occur significantly on the methyl 
group of methanol and not on its hydroxyl group, the 
[CH2DOH]/[CH30D] ratio in enhanced. The observed 
ratio of 10-20 can also be predicted but at a larger time, 
5 — 10 X 10** yr. We conclude that this class of models suc- 
ceeds to reproduce all the observations simultaneously, 
with the exception of the [CH2DOH]/[CH30D] ratio, for 
a density of about 5 x 10^ cm~'^ and at a time of 5000 

yr- 

4. DISCUSSION AND CONCLUSIONS 

Although multiply deuterated molecules in low mass 
protostars have been discovered and observed for a 
decade, published models have had difficulties in re- 
producing the observed abun dances, especially those of 
formaldehyde and methanol (jParise eTaDHool). Both 
species are believed to be mainly synthesised and to be 
observed on the grain surfaces during the cold and dense 
pre-coUapse phase (see lOberg et al.ir2011l ) and that they 
are observed in the gas when they sublimate off the 
grain mantles upon heating from the central star. On 
the grains, H2CO and CH3OH are thought to be the re- 
sult of the hydrogenation of iced CO. Their deuteration. 
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Fig. 3. — Mean deuterium fractionation of methanol and 
formaldehyde (solid lines) with their 1 sigma standard deviation 
(color levels) with time, obtained considering addition reactions 
only (§3.2). The h atched zones give the range of observed values 
dParise et al.ll200g'l . 




10' 10" 
Time (years) 

Fig. 4. — Mean deuterium fractionation of methanol and 
formaldehyde (solid lines) with their 1 sigma standard deviation 
(color levels) with time, obtained adding the abstraction and sub- 
stitution reactions (§3.3). The hatched zones give the range of 
observed values (Parise ct al. 2006 ). 
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therefore, depends on when exactly the two species are 
formed and how. We have shown that, based on our grain 
surface model GRAINOBLE, formation of H2CO and 
CH3OH from addition reactions alone fails to predict the 
observed deuterium fractionation. In contrast, if D and 
H abstraction and substitution reactions are added, the 
GRAINOBLE model can reproduce simultaneously the 
observed values. Therefore, these processes are crucial 
and more laboratory experiments and theoretical compu- 
tations are needed to better constrain their rates on the 
ices. Reproducing the observed [CH2DOII]/[CH3 0D] ra- 
tio remains a challenge, as it is still underestimated by 
the model. Previous studies have suggested that this 
may be caused by D an d H exchanges on the ic es during 
the sublimation phase (jRataiczak et al.' '2009') or acti- 
vated by photo lysis proc esses (Weber et al. 2009) or on 
the gas phase (Osamu ra et al.i 12004). all processes that 
wou ld be inefficient in alt ering the other isotopologues 
(see lRataiczak et al]|2011l for a detailed discussion). 

With the above exception, our model predicts the ob- 
served abundance ratios for high densities 5 x 10^ 
cm""^) and for relatively short times 5000 yr). We 
emphasise that this time corresponds only to the final 
stage at high density and not to the age of the conden- 
sation, which can be considerably larger. In fact, it is 



possible and even likely that the pre-s tellar cores spend 
a lon g time in a less dense phase (e.g. I Herein fc Tafallal 
l2007f) . However, the comparison of the observed II2CO 
and CH3OII deuteration with our model predictions sug- 
gests that the last phase at high density is short, just 
a few thousands years. In other words, soon after the 
central density reaches 5 x 10^ cm~'^ the protostellar 
collapse starts. This is a nice illustration of how chem- 
istry can help us to understand the past history of the 
protostar. 

In conclusion, our work leads to two important results: 
1) The pre-coUapse phase may last hundred thousands 
of years: however, when the central density reaches 
~ 5 X 10^ cm~'^, the collapse starts in about a few 
thousands years. 2) D and H atoms abstraction and 
substitution reactions are crucial in the grain surface 
chemistry and should be incorporated into models. 
Thus, more experimental and theoretical works are 
therefore needed to better contrain their efficiency 
and therefore the timescale needed to reproduce the 
observations. 



This work has been supported by I'Agence Nationale 
pour la Recherche (ANR), France (project FORCOMS, 
contracts ANR-08-BLAN-022). 
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